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ABSTRACT 

Context. New data from the Herschel Space Observatory are broadening our understanding of the physics and evolution of the outer 
regions of protoplanetary disks in star forming regions. In particular they prove to be useful to identify transitional disk candidates. 
Aims. The goals of this work are to complement the detections of disks and the identification of transitional disk candidates in the 
Lupus clouds with data from the Herschel Gould Belt Survey. 

Methods. We extracted photometry at 70, 100, 160, 250, 350 and 500 pm of all spectroscopically confirmed Class II members 
previously identified in the Lupus regions and analyzed their updated spectral energy distributions. 

Results. We have detected 34 young disks in Lupus in at least one Herschel band, from an initial sample of 123 known members in 
the observed fields. Using the criteria defined in Ribas et al. (2013) we have identified five transitional disk candidates in the region. 
Three of them are new to the literature. Their PACS-70 pm fluxes are systematically higher than those of normal T Tauri stars in the 
same associations, as already found in T Cha and in the transitional disks in the Chamaeleon molecular cloud. 

Conclusions. Herschel efficiently complements mid-infrared surveys for identifying transitional disk candidates and confirms that 
these objects seem to have substantially different outer disks than the T Tauri stars in the same molecular clouds. 

Key words, stars: pre-main sequence - protoplanetary disks - (stars:) planetary systems 


1. Introduction 


Transitional disks link the studies of planet formation and disk 
evolution around young stars. They are protoplanetary disks 
around young stars, optically thick and gas rich, with astronom¬ 
ical unit-scale inner disk clearings or radial gaps evident in their 
Spectral Energy Distributions (hereafter, SED). They lack excess 
emission at short-mid infrared (typically around 8- 12 pm ) but 
present considerable emission at longer wavelengths dMuzerolle 
|et al.||2010| ). Since their discovery by |Strom et al.|(|1989|) w ith 
IRAS data, the Spitzer Space Telescope ( [Werner et al.|2004] ) al¬ 
lowed the identification of a much larger population of these ob¬ 


jects ( |Cieza et al.|2010[ Merm et al.|2010[|Espaillat et al.|2014 ). 
It was possible to use the mid-infrared data to make good esti¬ 
mates of the dust distribution in the inner few astronomical units 
of the disks around T Tauri stars in nearby star-forming regions 
(e.g. Calvet et al.|2005[|Lada et al.|2006|). 


The Herschel Space Observatory ( [Pilbratt et al.|2010| grants 
us access to the outer disk region of these objects. Recent stud¬ 
ies suggest that the outer disks, as detected in far-infrared emis¬ 
sion, u ndergo substantial transformations in the transitional disk 
phase ( Cieza et al.|2011||Ribas et al.|2013| ). However, Herschel 
data are still lacking for most known transitional disks. 

The main objective of this article is to revisit and extend our 
knowledge of known Class II objects in the Lupus dark clouds, 


* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 


and to identify transitional disk candidates in the region using 
Herschel photometry, complemented with photometry from pre¬ 
vious studies from the optical to the mid-infrared. With Herschel 
we have the opportunity to study objects with large inner holes, 
reaching for the furthest regions of their disks, where Spitzer was 
unable to access with the required sensitivity (IRAC and IRS had 
different wavelength ranges, between 3.5 and 8.0 pm and 5 and 
35 pm, respectively, and the sensitivity of MIPS, with pass bands 
between 24 and 160 pm, was limited). Therefore, we are not only 
able to study the outer regions of known transitional disks, but 
also detect new candidates, with larger cavities and outer disks. 

This work is organized as follows: §[2]describes the Herschel 
observations and data reduction and the procedure applied to 
identify and extract the fluxes for the known objects. In §[3]we 
identify new transitional disk candidates, describe the procedure 
to build the SEDs and discuss some objects with certain peculiar¬ 
ities. In § [4] we discuss the properties of the newly found transi¬ 
tional disk candidates and in § [5] we present the conclusions of 
this work. 


2. Observations and data reduction 


The Lupus dark clouds were observed by Herschel as part of 
the Herschel Gould Belt Survey (HGBS hereafter, 


2010). It is one of the star-forming regions located within 150 


Andre et al. 


to 200 pc from the Sun and has a mean age of a few Myrs (see 
Comer6n||2008| for a review about the region). It contains four 


main regions, with one of them, Lupus III, being a rich T Tauri 
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association. Located in the Scorpius-Centaurus OB association, 
the massive stars in the region are likely to have played a signif¬ 
icant role in its evolution._ 

By employing the PACS ( [Poglitsch et al.||2010| ) and SPIRE 
( Griffin et al.|2010| ) photometers in parallel mode, Herschel has 
the capability for large scale mapping of star forming regions 
in the far-infrared range. Early Herschel observati ons o f Lupus 
from the HGBS were already published by |Rygl et al.| ( [2013] >, 
where the observing strategy and main results are described. 

We produced our maps using Scanamorphos (Roussel |2012| ) 
version 21.0 for the PACS maps and version 18.0 for the SPIRE 
ones. The obsids of the Lupus Herschel parallel mode obser¬ 
vations processed for this work are (1342-) 189880, 189881, 
203070, 2203071, 2203087, 2203088, 2213182 and 2213183 
from program KPGT_pandre_l. The observations were made in 
parallel mode for 70, 160,^50, 350 and 500 pm, and in prime 
mode for 100 pm ( |Rygl et al. 12013] ). This is the reason why some 
of the sources do not have images at 100 pm in Figures [5] to [9j 
A complete census of the sources detected with Herschel will 
be published by the Herschel Gould Belt Survey consortium in 
an upcoming paper (Benedettini et al. in prep.). 


2.1. Point source photometry 


Our first objective was to detect and classify sources in the im¬ 
ages obtained with PACS and SPIRE. We used coordinates of 
spectroscopically confirmed Class II members from different 
sources, compilinga list of 2 17 targets. This compilation was de¬ 
veloped by |Ribas et al. ( [2014 L where these objects were selected 
for having known extinction values and spectral types from op¬ 
tical spectroscopy. Of this sample of 217 objects, only 123 are 
found in the fields observed by Herschel. 

We applied the sourceExtractorSussextractor algorithm in 
our maps, with the Herschel Interactive Processing Environment 
(HIPE), version 12.0, using a detection threshold of S/N> 3. 
Then we cross-matched the resulting list of detected objects with 
our pre-selected targets, using a search radius between 5” and 
10”. We selected only those objects with detections in PACS. 
After this process, we found 37 objects with at least one de¬ 
tection in one PACS band, 35 in the 70 pm maps and 2 in the 
100 pm maps. We then validated all of them by visual inspec¬ 
tion. This way we were able to double check our detections. The 
background emission becomes more significant at longer wave¬ 
lengths, where detections are in general less frequent, or false 
detections appear, due to bright filaments and ridges. This is the 
reason why we selected only PACS-detected objects. The com¬ 
bination of quantitative and qualitative analysis was essential to 
produce reliable detections and photometry. After this process 
we discarded 4 objects as not having clear detections in the maps 
and added 2 objects that were found by visual inspection. They 
are flagged in Table [5] 

After this process, 35 objects of the list had one clean detec¬ 
tion in at least one Herschel PACS or SPIRE band. 

We applied the sourceExtractorDaophot algorithm to ex¬ 
tract the photometry for PACS, and the annular SkyAperturePho- 
tometry task for SPIRE. All point-source fluxes were aperture 
corrected with the recommended aperture corrections for both 
instruments (see the PACS Photometer Point Spread Function 
Technical Note from April 2012, and Sect. 6.7.1 of the SPIRE 
Data Reduction Guide of HIPE 12.0). Table [T] shows the set of 
aperture, inner and outer sky radii for each Herschel band, as 
well as the corresponding aperture corrections. 

Table [2] shows the sensitivity limit for each Herschel band, 
as stated in the SPIRE/PACS Parallel Mode Observers’ Manual 


Table 1: Aperture, inner and outer radii used in the photometry 
extraction process for each band, with their corresponding aper¬ 
ture corrections. 


Band 

R aperture ( ) 

R inner ( ) 

R outer ( ) 

Correction 

PACS 70 

6 

25 

35 

1.5711 

PACS 100 

6 

25 

35 

1.6804 

PACS 160 

12 

25 

35 

1.4850 

SPIRE250 

22 

60 

90 

1.2796 

SPIRE350 

30 

60 

90 

1.2396 

SPIRE500 

42 

60 

90 

1.2937 


(section 2.3). The values correspond to a nominal scan direction 
and a scan speed of 60”/s. 


Table 2: Sensitivity limits for each Herschel band in paral¬ 
lel mode at 60”/s, as stated in the PACS/SPIRE Parallel Mode 
Observers’ Manual. 


Filter 

Sensitivity (mJy/beam) 

PACS 70 

21.0 

PACS 100 

24.7 

PACS 160 

47.0 

SPIRE 250 

12.6 

SPIRE 350 

10.5 

SPIRE 500 

15.0 


Upper limits for non-detections were derived by adding 
scaled synthetic Point Spread Functions (PSFs) at the expected 
position of the sources until the sussExtractor detection algo¬ 
rithm did not recover them. In most cases visual inspection 
was applied. The calibration errors for PACS and SPIRE are 
5% and 7% (as stated in their observer manuals), respectively. 
To keep consistency with Ribas et al. (2013) we chose a more 
conservative estimation, 20% for SPIRE as error values, taking 
into account that the background emission becomes increasingly 
stronger at longer wavelengths. In the case of PACS, after com¬ 
paring several photometric extraction algorithms ( sussExtractor , 
Daophot , a nnularSkyAperture Photometry and Hyper (Traficante 
et al.|2014| ), and extraction apertures with the MIPS70 fluxes of 
a subsample of clean sources, we selected a conservatively er¬ 
ror value of 25% of the flux, to account for the different flux 
values obtained. Table [5] gives the stellar parameters of these de¬ 
tected sources, as given in the literature, and Table [6] presents 
their Herschel point source fluxes and upper limits. As stated 
before, at least one point source was detected for at least one 
wavelength of Herschels instruments for all objects in Table [6] 
(see Figures 5 to 9). 

We then looked for bright Spitzer 24 pm and 70 pm sources 
around the targets, to check for possible contamination from 
other nearby mid-infrared sources. We used a radius of 42”, 
as this was the value we used for the largest wavelength band 
when extracting the Herschel photometry (which corresponds to 
SPIRE 500 pm). The results of this analysis are shown in Tables 
[3] and [4] 

In the cases where MIPS70 sources were found near our 
sources we inspected the PACS70 images of the targets. Except 
for sources Sz65, Sz66, Szl03, Szl04 and 2MASSJ1608537 (as 
explained below), we did not found evidence of the presence of 
extra sources inside the matching radius (see Figures [5] to [9]). 
This means that Herschel did not detect them, so we assume 


2 
































I. Bustamante et al.: Identification of new transitional disk candidates in Lupus with Herschel 


that the full contribution to the flux measured at PACS70 comes 
from our target sources, and not from contaminants. In Table]?] 
we maintain all the sources with possible contamination for ref¬ 
erence. 

Objects Sz65 and Sz66 form a pair, as is explained in Section 
|3.3.6 and this is taken into account in their photometry extrac¬ 
tion. The MIPS excesses seen for objects Szl03 are related with 
object Szl04, and vice versa. We only flagged objects with po¬ 
tential contaminations in case they were detected (see Table [6]). 
When they are not detected and upper limits are calculated in¬ 
stead, we do not take into account the possible contaminations. 

The nearby source near object 2MASSJ1608537 presented 
significantly larger MIPS24 and MIPS70 fluxes than our target. 
In this case, we concluded that the infrared excesses measured 
in the Herschel image corresponded to this other object instead 
of our target, and thus discarded object 2MASSJ1608537, as 
its Herschel photometry could not be unequivocally associated 
with it. The other object is a known YSO candidate, SSTc2d 
J160853.2-391440. It is probably another member of the Lupus 
region, yet not confirmed spectroscopically and therefore not in¬ 
cluded in our study. 

In other cases other contaminant sources were found, but 
with sufficiently weaker MIPS24 and MIPS70 fluxes as not to 
affect our measurements. These sources are marked in Tabled 
In these cases we assumed that the Herschel fluxes measured 
correspond to our targets as the major contributors. 

Another interesting object is Szl08B, which presents two 
possible contaminants that could affect its photometry. The one 
furthest away, at 38.8”, has larger MIPS24 flux. The closest one, 
at 9.7”, has a higher MIPS70 flux. Thus, the Herschel photom¬ 
etry of Szl08B could be contaminated by the closest one and 
therefore it is flagged as unreliable in Table]?] 

After this analysis, our final set of sources decreased to 34. 
We maintain the image of object 2MASSJ1608537 in Figure[7] 
for completeness. 


3. Results 

3.1. Identification of transitional disk candidates 


To identify transitional disk candidates we applied the procedure 
developed in [Ribas et al. (2013), which defines a transitional disk 
candidate as having a K -u < 0 and a 12-70 > 0, where ct\ 1 -\ 2 is 

defined as a kl -x 2 = — log^n-iog^) > and A 1S measured in 


and in erg • s _1 • cm -2 . 

The K-band flux comes from 2MASS ( [Skrutskie et al.|2006| ), 
the 12 jum flux from WISE and the 70 ji m one from PACS. 
Some exceptions to this rule were applied when these fluxes 
were unavailable for certain objects. Sz95 has no flux at 70 yum, 
so we used PACS 100 yum instead. We did the same for 2MASS 
J1608149. An exceptional case i s obje ct SSTc2dJ160111.6- 
413730, discussed later on section [3?3.5| as no PACS flux was 
available, and therefore could not be classified. 

With this procedure we detected five objects that satisfy the 
transitional disk candidate criteria: SSTc2dJ161029.6-392215, 
Sz91, Szlll, Szl23 and 2MASS J1608149. In Figure^we rep¬ 
resent the objects that fulfilled our selection criteria. Objects 
Sz91 and Szlll were previously known transitional disks from 
Spitzer data ( |Merm et al.|2008| ), and they occupy a different re¬ 
gion in the diagram than the others. They present little to no 
excess at 12 yum and a steep rise between 12 and 70 yum in their 
SEDs (see Figure [2]). This could be indicative of a large cavity 
in both cases, mostly empty of small dust grains. In the case of 


Sz91, a resolved image of this inner hole has already been ob¬ 
tained in the sub-millimeter with SMA ( Tsukagoshi et al.|2014 ), 
being one of the largest imaged so far (~65 AU) in a T Tauri 
disk. Szlll is likely to have a similarly large and empty cavity. 


The rest of the objects present excesses at intermediate wave¬ 
lengths, as well as fulfilling our color criteria (see Figure [2] for 
a representation of their SEDs). Objects 2MASS J1608149 is a 
special cases, as explained at subsection |3. 3.3 

In Figure [I] three additional objects present error bars falling 
in the transitional disk detection criteria: SSTc2dJ160703.9- 
391112, RXJ1608.6-3922 and Szl08B. Their SEDs present far 
infrared excesses, but mid-infrared ones as well. As they do not 
strictly fulfill our selection criteria, they are not to be consid¬ 
ered transitional disk candidates in this work, yet they deserve 
attention. 


Analysis on the remaining sample of T Tauris in the Lupus 
regions is outside the scope of this wo rk, a lthough some special 
cases are discussed in some detail in § |3.3| 


3.2. Spectral Energy Distributions 


To construct the SEDs we used ancillary data from g round -based 
optical, near-infrared, Sp itzer and WISE data (Mortier et al 


|2011[ |Comer 6n 2008; [Allen et aL]|2007[ |Krautter et al.||1997] T 
These photometric data provide continuous wavelength cover¬ 
age between 0.4 and 500 yum. Figures [2] [3] and]?] show the SEDs 
of the transitional disk candidates and the rest of the detected 
objects, respectively. We used the interstellar extinction law of 
Weingart ner & Draine ( 200l] ). The A v values were extracted 


from Mortier et al.| ( [2011| ) and |Merin et aL| ( |2008| ), except in 
7 sources, identified with an asterisk in Table [6] for which we 
derived them from the observed optical and near-infrared pho¬ 
tometry and their spectral types. 


We inspected all the images in 2MASS band J of the entire 
sample of detected objects. All of them were confirmed as point 
sources except those marked in Table [5] with a j , where minor 
extended emission cannot be discarded. 


Fig. 1: Identificatioji of transitional disks candidates applying the 
color criteria from Ribas et al. ( [2013 ) to the Lupus Herschel data. 
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Table 3: Targets with possible contamination within the used apertures in MIPS24. The second column shows the fraction of total 
MIPS24 flux in the aperture with that of the target. 


Object 

Esource/Eaper 

N 

Companion distance (”) 

Band 

Sz65' 

0.4 

1 

6.5 

PACS70 

Sz66’ 

2.9 

1 

6.2 

PACS70 

SzlOO 3 

0.2 

1 

23.4 

SPIRE350 

Szl03 2 

0.8 

10 

41.2 

SPIRE500 

Szl04 2 

2.7 

10 

36.4 

SPIRE500 

Szl08B 

1.4 

3 

38.8 

SPIRE500 

2MASSJ1608537* 

9.3 

1 

6.9 

PACS70 

Szll2 

0.5 

4 

31.9 

SPIRE500 

SSTc2dJ161029.4-392215 3 

0.1 

2 

7.4 

PACS70 

SSTc2dJ160002.4-422216 3 

0.4 

2 

39.8 

SPIRE500 

SSTc2dJ160111.6-413730 3 

0.2 

1 

24.2 

SPIRE350 


Notes. 1 The MIPS24 contamination affecting Sz65 is caused by Sz66, and vice versa, as can be seen in Figure [ 5 ] This is taken into account in 
the photometry analysis. 

2 The MIPS24 contamination affecting Szl03 is caused by Szl04, and vice versa, as can be seen in Figures [6] and [ 7 ] This is taken into account in 
the photometry analysis. 

3 Objects with sufficiently low contamination from nearby sources as to not affect our Herschel measurements. 

* This is the only object that has been removed from the list of detected sources, because its Herschel photometry cannot be unequivocally 
associated to it. 


Table 4: Same as Table[3] but with MIPS70 possible contamination. 


Object 

Esource/Eaper 

N 

Companion distance (”) 

Band 

Sz66 

A6%. mJy 

1 

6.2 

PACS70 

Sz98 2 

0.2 

1 

26.3 

SPIRE350 

SzlOO 2 

0.7 

1 

32.3 

SPIRE500 

Szl03> 

1.4 

1 

23.2 

SPIRE350 

Szl04> 

0.7 

1 

23.2 

SPIRE350 

Szl08B 

212. mJy 

1 

9.7 

PACS70 

2MASSJ1608537* 

U9.0 mJy 

1 

6.9 

PACS70 

Szll2 2 

0.7 

1 

28.6 

SPIRE350 


Notes. 1 The MIPS70 contamination in Szl03 is caused by Szl04, and vice versa, as can be seen in Figure [6] and [ 7 ] This is taken into account in 
the photometry analysis. 

2 Objects with sufficiently low contamination from nearby sources as to not affect our Herschel measurements. 

* This is the only object that has been removed from the list of detected sources, because its Herschel photometry cannot be unequivocally 
associated to it. 

mJy In these cases the MIPS photometry of the target was not available, but the one of the companions was. Therefore, we present the actual flux 
of the objects inside the search radius, instead of the percentage. 


3.3. Individual objects of interest 
3.3.1. Sz91 and Szlll 


Sz91 and Szlll had been previously identified as potential tran 
sitional disks based on Spitzer data in [Merin et al.| ( f2008| and 


have large H a equivalent widths of 95.9 and 145.2 A (Hughes 
et al.|1994j ), which indicate youth and active gas accretion to the 


stars. 

The optical images of S z91 reveal a known companion at 
approximately 8" ( Ghez et al. 1997), which corresponds to 1200 
to 1600 AU at probably the distance of the star. However, the 
centroids of all the Herschel -detected sources match the coordi¬ 
nates of the primary source, Sz91 (see Figure |5j, so we assume 
no contribution from the companion to the measured fluxes. 


stead of the wide ones found at Sz91 and Szlll. They also 
present small excesses at near-infrared wavelengths (3-10 pm), 
that could signal the presence of an inner optically thick disk. 

3.3.3. 2MASS J1608149 

Since object 2MASS J1608149 was not detected on the PACS70 
band, we then used the PACS100 one to classify it. The target 
then fulfilled the transitional disk detection criteria. As its SED 
shows [2j the Herschel far-infrared excess is not as significant as 
to definitely classify the object as a transitional disk. 

A further analysis including modeling could reveal the true 
nature of this object. As for this work, we maintain object 
2MASS J1608149 as a transitional disk candidate. 


3.3.2. SSTc2dJ161029.6-392215 and Sz123 

Objects SSTc2dJ161029.6-392215 and Szl23 present SEDs 
characteristic of transitional disks but with narrow gaps, in- 


3.3.4. Sz68 

Figure [5] shows an extended feature at all wavelengths to the 
North of the Sz68 object. This is an already known Herbig-Haro 
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object driven by Sz68 ( Moreno-Corral et al.|1995| ). Interestingly, 
its emission dominates at wavelengths longer than 100 pm, 
which probably contaminates the SED of this object with un¬ 
related flux to the star itself. A similar effect was discovered 
around the star T54 in Chamaeleon I ( [Matra et al.l|2012| ), er¬ 
roneously classified before as a transitional disk due to nearby 
extended emission (see also |Lestrade et al.|2012[ ). 


3.3.5. SSTc2dJ160111.6-413730 

SSTc2dJ160111.6-413730 has a very peculiar SED (see Fig. [4]) 
with no excess at any wavelengths between optical and 160 pm 
but substantial potential excesses at 250 to 500 /mi associated to 
a clear point source in the images (see Fig. [9]). The low 70 /mi 
upper limit is nevertheless consistent with the sensitivity limit of 
the observations. 

Previous claims of such extreme type of objects have been 
made in the context of the DUNES Herschel Key Programme, 
who named these obj ects ’’cold disks” t o make reference to 
their extreme nature ( [Krivov et al.||2013] ). Other authors have 
pointed towards background galaxies as the possible explana¬ 
tion ( [Gaspar & Rieke||2014| ). It did not fulfill our classification 
criteria as transitional disk candidate. 


all spectroscopically confirmed Class II members of the associ¬ 
ation. Two transitional disks were known from previous works 
(Sz91 and Szl 11) but the other three are new to the literature. 

The Spitzer study of the Lupus clouds found a global disk 
fraction of 50 - 60% ( [Ribas et al.||2014| ), compatible with a 


relatively young age for the region of 1 - 2 Myrs (Comeron 
2008[ ). However, the total sample in |Merfn et al.| ( [2008| ) did in 


elude objects without spectroscopic membership confirmation 
so it could have certain level of contaminatio n from background 
galaxies and/or old dusty AGB stars (see Mortier et al. 2011 [ for 
a spectroscopic survey of Spitzer selected candidate members 
in Lupus). This work however uses an input target list where 
all objects have spectroscopic membership confirmation and is 
therefore not affected by that. 

The observed fraction of transitional disks in Lupus based 
on our Herschel detected sample is -15%, which is compara¬ 
ble to the fractions measured in other such young and nearby 
regions |Espaillat et al . 2014| ). This corrects a surprisingly low 
incidence of this type of objects from the Spitzer sample re¬ 
ported in |Merfn et al. ( 2008| , where only 2 objects were iden¬ 
tified as transitional, from a photometric sample of 139 Class 
II and III sources. This demonstrates how Herschel efficiently 
complements mid-infrared surveys for the specific study of tran¬ 
sitional disks. 


3.3.6. Sz65 and Sz66 

Objects Sz65 and Sz66 form a binary system, with a separa¬ 
tion of 6.4” (Lommen et al. |2010 ), unresolved by Herschel (see 
Figure [5]). We inspected their SEDs to select which one of them 
was the main contributor. We deduce the primary source to be 
Sz65 by studying the contour diagram of the image. The con¬ 
tour lines were centered at the position of source Sz65. Their 
MIPS24 and MIPS70 fluxes support this association. Then we 
set the Herschel fluxes of Sz66 as upper limits in Tableland in 
Figure [3] 


4. Discussion 


4.1. Detection statistics 

From an initial sample of 217 known Class II sources with spec¬ 
troscopic membership confirmation in the Lupus clouds, 123 ob¬ 
jects fell in at least one of the fields observed by the Herschel 
Gould Belt Survey. In particular, 92 and 98 were in the PACS 
and SPIRE Lupus III maps, respectively; 12 and 15 were found 
in Lupus I PACS and SPIRE mosaics; and 10 objects were ob¬ 
served in Lupus IV with both instruments. From these sources, 
we detected 35 objects, 32 in PACS70, 2 in PACS 100 and 1 in 
SPIRE. We then discarded one, 2MASSJ1608537, as we could 
not unequivocally assign its Herschel photometry to the target. 

Compared with a similar study in Chamaeleon ( Ribas et ah] 
2013| ), here we report a Herschel detection rate of - 28% of 


the previously known objects in the region, comparable with the 
percentage in that work (-30%). The small different distances 
to the regions (~ 150 pc for Chamaeleon vs - 200 pc to Lupus 
III) might account for the small difference in the detection rate. 
Further analysis of the detection statistics of known and new 
sources in these regions with the HGBS data will be presented 
by Benedittini et al. (in prep.). 


4.2. Incidence of transitional disks in Lupus 

We report the detection of five transitional disk candidates in the 
Lupus molecular clouds based on Herschel and previous data of 


4.3. Brighter PACS-70 fluxes in transitional disks 


Recent studies have reported brighter 70 pm fluxes in transi¬ 
tional disks as compared with the median S ED of ’’standard” 
T Tau ri stars in the same regions ( [Cieza et al.|2011| |Ribas et aL 


2013). To check whether this phenomenon is dependent on the 


conditions of certain molecular clouds, we compared the SEDs 
of the detected disks in Lupus with the median SED of T Tauri 
stars in the region. Figures [2] [3] and [4] show the median SED nor¬ 
malized to the dereddened J-band flux for comparison, and Table 
[ 7 ] shows the fluxes of this median SED. 

As shown in Figure [2] most transitional disk candidates iden¬ 
tified in section show the same phenomenon: the PACS 
fluxes are systematically above the median SED of the T Tauri 
stars. Therefore, the result is confirmed so far around the isolated 
transitional disk T Cha ( [Cieza et al.|2011| ), in the Chamaeleon I 
cloud ( [Ribas et al.|2013| ) and in the Lupus molecular clouds (this 
work). 

The interpretation of this phenomenon requires detailed SED 
modeling with radiative transfer disk models. However, the 
higher PACS fluxes as compared with the median SED of the 
T Tauri stars suggests that the inner and outer disks of these ob¬ 
jects follow different evolutionary paths by the time the inner 
gaps or holes are formed in the inner disks. Thus, the evolution 
of both inner and outer disk regions might not be dynamically 
decoupled. 


4.4. SED population analysis 

Merm et aL| ( [2008| classified the YSOs population in Lupus by 
studying the shape of their SEDs, and comparing the median 
SED of the CTTS from Taurus with their data. They grouped 
the objects in four categories, based on whether their slope de¬ 
cay like a classical accreting optically thick disk around a low- 
mass star (T-Type), present infrared excess clearly weaker (L- 
Type) or stronger (H-Type) than the median SED or where no 
excess is detected (E-Type). They also introduced the category 
cold disks (LU-Type), which corresponds to our criteria for tran- 
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sitional disks. In their work, they made use of Spitzer and ground 
based photometry. 

We did the same analysis for the 34 objects detected by 
Herschel , but now also taking into account the photometry we 
extracted. For simplification, we used a different nomenclature, 
grouping our objects in only three categories to unveil their evo¬ 
lutionary state: higher infrared excess than the median values 
(H), indicating an early evolutionary stage; weaker (W), showing 
a more evolved one; and transitional candidate (T). The results 
can be found in Table [5j where we also comp are our classifica¬ 
tion to the aforementioned work by |Merfn et al. (2008). Another 
difference is that instead of using the Taurus median SED, we 
have used the Lupus one. 

Apart from the transitional nature for the three new objects 
found in this work (~ 10%), we find that by adding HerscheVs 
photometry to the SEDs we can update the type found in |Merfn| 
[et al.||200S] ) for another 9 objects (~ 27%). The final population 
for our sample is then classified as follows: ~ 56% primordial 
disks (this groups T-Tauri types and stronger infrared excesses), 
~ 29% of evolved disks and ~ 15% for transitional disk candi¬ 
dates. This compares well to the 67%, 26% and 7% fractions 
for the same types of objects from the Spitzer- only study. Table 
[5] shows the Spitzer- only and Spitzer+Herschel classifications. 
The relative overall fractions have not changed substantially, ex¬ 
cept on the larger fraction of transitional disk candidates. 
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5. Conclusions 


We have detected 34 objects in the Lupus cloud using Herschel 
and broadened the study of the YSOs in the association. Having 
confirmed that our Herschel data enable improved characteriza¬ 
tion of outer regions of protoplanetary disks, we have identified 
additional transitional disks in our sample. The fraction of tran¬ 
sitional disks detected in our study, comparable in size to Ribas 
et al. ( 2013| ), corrects the relatively low one from Spitzer de¬ 
tections. Also, we find that all these objects have 70 pm fluxes 
brighter than the median SED, showing perhaps another possible 
intrinsic property of them. Finally, we have carried out a popula¬ 
tion analysis by studying the SED shapes of the detections, and 
updated the morphological type of several sources. Further stud¬ 
ies and modeling will give us more information of these objects. 
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Table 5: Class II objects detected by Herschel in the Lupus clouds. The object types stand as higher than the median infrared excess 
(H), weaker (W) and transitional disk candidates (T). See text for details. 


Object 

R-A.72000 

Dec/2000 

A v (mag) 

SpT 

Type (, Spitzer **) 

Type (This Work) 

References'" 

Lupus I 

Sz65 

15:39:27.8 

-34:46:17.1 

1.0 

M0 

H 

H 

2 

Sz66 v 

15:39:28.3 

-34:46:18.0 

2.0 

M3 

H 

H 

2 

Sz68 

15:45:12.9 

-34:17:30.6 

1.5 

K2 

H 

H 

2 

Sz69 

15:45:17.4 

-34:18:28.3 

4.0 

Ml 

W 

W 

2 

Lupus III 

SSTc2dJ160703.9-391112 

16:07:03.8 

-39:11:11.3 

2.2 

M5.5 

- 

H 

1 

Sz90 

16:07:10.1 

-39:11:03.3 

5.0* 

K8 

W 

H 

2,5 

Sz91 

16:07:11.6 

-39:03:47.5 

1.0 

M0.5 

T 

T 

2 

Sz95 

16:07:52.3 

-38:58:05.9 

1.0 

M1.5 

W 

W 

2 

Sz96 

16:08:12.6 

-39:08:33.5 

1.4 

M2 

H 

W 

1 

2MASSJ1608149 

16:08:14.9 

-38:57:14.6 

7.0* 

M4.7 

H 

T 

3,5 

Sz98 

16:08:22.5 

-39:04:46.5 

2.5 

K5 

H 

H 

1 

SzlOO 

16:08:25.8 

-39:06:01.2 

1.2 

M4.5 

H 

H 

1 

Szl02 

16:08:29.7 

-39:03:11.0 

2.5 

K2 

- 

H 

1 

Szl03 

16:08:30.3 

-39:06:11.2 

1.0 

M4.5 

H 

H 

1 

Szl04 

16:08:30.8 

-39:05:48.9 

0.7 

M5.5 

H 

H 

1 

HR5999 

16:08:34.3 

-39:06:18.2 

1 .0* 

A7 

H 

W 

2,5 

RXJ1608.6-3922 

16:08:36.2 

-39:23:02.5 

3.0* 

K6 

- 

H 

4,5 

Szl08B 

16:08:42.9 

-39:06:14.4 

0.8 

M5.5 

W 

W 

1 

SzllO 

16:08:51.6 

-39:03:17.7 

0.2 

M3 

H 

H 

1 

Szlll 

16:08:54.6 

-39:37:43.1 

0.0 

M1.5 

T 

T 

2 

Szll2 7 

16:08:55.5 

-39:02:33.9 

1.0 

M4 

W 

W 

2 

SSTc2dJ160901.4-392512 

16:09:01.4 

-39:25:11.9 

0.6 

M3.5 

H 

H 

1 

Szll4 7 

16:09:01.9 

-39:05:12.4 

1.3 

M4 

H 

H 

1 

Szll7 

16:09:44.4 

-39:13:30.1 

1.0 

M2 

H 

W 

2 

Szll8 7 

16:09:48.7 

-39:11:16.9 

2.6 

K7 

H 

H 

1 

SSTc2dJ161029.6-392215 

16:10:29.6 

-39:22:14.5 

0.0 

M5.5 

H 

T 

1 

Szl23 

16:10:51.6 

-38:53:13.8 

0.0 

M3 

W 

T 

2 

Lupus IV 

SSTc2dJ160002.4-422216 

16:00:02.4 

-42:22:14.6 

2 .0* 

M3.5 

H 

H 

1,5 

IRAS 15567-4141 J 

16:00:07.4 

-41:49:48.4 

2.0 

M6.5 

W 

W 

1 

Szl30 J 

16:00:31.1 

-41:43:36.9 

0.60 

M2 

H 

W 

1 

F403 7 

16:00:44.5 

-41:55:31.0 

2.0 

K0 

W 

H 

2 

Szl3L 

16:00:49.4 

-41:30:03.9 

0.0 

M2 

H 

H 

2 

SSTc2dJ160111.6-413730 v 

16:01:11.5 

-41:37:29.9 

0.0 

M8.5 

W 

W 

1 

Szl33 7 

16:03:29.4 

-41:40:01.8 

10 .0* 

K2 

H 

H 

1 


Mortier et al. 


( 2011 ) , ( 2 ) 


Comeron 


42008) , (3) 


Allen et al. 


(2007 ), (4) 


Krautter et al. 


(1997) 


Notes. (r) References for the spectral type and A y : (1) 
and (5) this work. 

(i) The objects marked with j presented potential small extended emission at J band. Nevertheless, all of them have spectroscopic confirmation of 
being young stellar objects. 

(v) These objects were missed by our detection method and were detected in the visual inspection. 

( * } Extinction values denoted with an asterisk indicate that were computed in this work either because they were not available in the literature or 
the onesfound did not provide appropriate fits to the model photospheres. 

20081. 


(**) 


Merin et al. 
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Table 6: Herschel photometry for the 34 YSOs detected in the maps. When no source was detected at a certain band, an upper limit 
is provided. 


Object 

E70 /um (Jy) 

FlOO^m (Jy) 

El 60 fim (Jy) 

E250 um (Jy) 

F350^m (Jy) 

F 5 00^m (Jy) 

Lupus I 

Sz65 

0.47+0.12 

- 

0.55+0.14 a 

0.77+0.15 a 

0.95+0.19 a 

<1.59 

Sz66 

0.45+0. ll a 

- 

0.55+0.14 a 

0.81+0.16 a 

0.98+0.20 a 

<1.63 

Sz68 

3.01+0.75 a 

3.80+0.95 a 

11.77+2.94 a 

13.06+2.61 a 

10.22+2.04 a 

6.74+1.35 a 

Sz69 

0.12+0.03 

0.19+0.05 

<0.31 

<1.59 

<0.88 

<2.13 

Lupus III 

SSTc2dJ160703.9-391112 

0.07+0.02 a 

0 .10+0.02 

0.19+0.05 a 

<0.33 

<0.44 

<0.43 

Sz90 

0.36+0.09 

0.37+0.09 a 

0.26+0.07 a 

0.51+0.10 a 

0.50+0.10 a 

<0.61 

Sz91 

0.51+0.13 

0.68+0.17 a 

0.72+0.18 a 

0.86+0.17 a 

0.62+0.12 a 

0.38+0.08 a 

Sz95 

<0.05 

0.05+0.01 a 

0.09+0.02 

<0.28 

<0.30 

<0.25 

Sz96 

0.14+0.03 a 

0.12+0.03 

0 .10+0.02 

<1.49 

<1.87 

<1.96 

2MASSJ1608149 

<0.08 

0.07+0.02 a 

<0.11 

<0.48 

<0.25 

<0.21 

Sz98 

0.64+0.16 a 

0.58+0.14 a 

0.50+0.12 a 

0.93+0.19 a 

1.14+0.23 a 

<2.12 

SzlOO 

0.17+0.04 

0.23+0.06 a 

0.29+0.07 a 

0.59+0.12 a 

<1.99 

<1.33 

Szl02 

0.34+0.09 

0.37+0.09 a 

0.22+0.06 a 

<3.25 

<4.07 

<4.37 

Szl03 

0.11+0.03 a 

0.16+0.04 a 

<0.23 

<2.13 

<1.64 

<1.20 

Szl04 

0 .10+0.02 a 

0.18+0.04 a 

<0.19 

<2.08 

<1.97 

<2.13 

HR5999 

4.25+1.06 

3.20+0.80 

1.60+0.40 a 

0.92+0.18 a 

0.47+0.09 a 

<3.26 

RXJ1608.6-3922 

0.40+0.10 a 

0.69+0.17 a 

1.36+0.34 a 

2.55+0.51 a 

2.83+0.57 a 

2.40+0.48 a 

Szl08B 

0.20+0.05 a 

0.21+0.05 a 

<0.61 

<5.81 

<9.68 

<8.85 

SzllO 

0.13+0.03 a 

0.09+0.02 a 

<0.26 

<0.57 

<0.48 

<0.76 

Szlll 

1.16+0.29 a 

- 

1.73+0.43 

1.31+0.26 

0.94+0.19 

0.58+0.12 a 

Szll2 

0.10+0.03 a 

0.04+0.01 

<0.16 

<0.33 

<0.27 

<0.38 

SSTc2dJl 60901.4-392512 

0.08+0.02 a 

0.09+0.02 a 

0.12+0.03 a 

0.15+0.03 a 

0.13+0.03 a 

0.07+0.01 a 

Szll4 

0.28+0.07 

0.26+0.06 a 

0.13+0.03 a 

<0.53 

<0.51 

<1.11 

Szll7 

0.06+0.02 a 

0.11+0.03 a 

0.17+0.04 a 

<0.82 

<0.70 

<0.65 

Szll8 

0.35+0.09 

0.46+0.11 

0.37+0.09 

0.49+0.10 a 

0.56+0.11 a 

<0.98 

SSTc2dJ161029.6-392215 

0.07+0.02 a 

0 .10+0.02 a 

<0.13 

<0.29 

<0.27 

<0.26 

Szl23 

0.27+0.07 a 

0.37+0.09 a 

0.45+0.11 

0.17+0.03 a 

0 .12+0.02 a 

0.08+0.02 a 

Lupus IV 

SSTc2dJ160002.4-422216 

0 .10+0.02 a 

0.20+0.05 a 

0.21+0.05 a 

0.22+0.04 

0.23+0.05 a 

0.17+0.03 a 

IRAS15567-4141 

0.11+0.03 a 

0.05+0.01 a 

<0.11 

<0.38 

<0.34 

<0.22 

Szl30 

0.11+0.03 

0.12+0.03 a 

0.05+0.01 a 

<0.39 

<0.34 

<0.34 

F403 

1.04+0.26 a 

1.26+0.3E 

1.57+0.39 a 

1.33+0.27 a 

0.93+0.19 a 

0.58+0.12 a 

Szl31 

0.11+0.03 

- 

<0.10 

<0.24 

<0.35 

<0.25 

SSTc2dJ160111.6-413730 

<0.05 

- 

<0.09 

0 .11+0.02 a 

0.07+0.01 

<0.16 

Szl33 

0.21+0.05 a 

0.38+0.10 

0.32+0.08 

0.37+0.07 a 

0.30+0.06 a 

0.20+0.04 a 


Table 7: Normalized flux densities of the median SED, lower SED (first quartile) and upper SED (fourth quartile) of the Class II 
objects in Lupus I, III and IV. The number of detections in each band is also included. The transitional disk candidates are excluded 
from this count. All values are normalized to J band. 


Band(//m) 

Median 

First Quartile Fourth Quartile 
(F^ arbitrary units) 

Detections 

J 1.25 

1.000 

1.000 

1.000 

28 

H 1.62 

1.048 

0.982 

1.186 

28 

K 2.15 

0.761 

0.673 

1.015 

28 

IRAC 3.6 

0.331 

0.260 

0.501 

23 

IRAC 4.5 

0.259 

0.182 

0.385 

24 

IRAC 5.8 

0.193 

0.143 

0.284 

24 

IRAC 8.0 

0.167 

0.111 

0.223 

25 

MIPS 23.68 

0.090 

0.054 

0.122 

27 

PACS 70 

0.049 

0.025 

0.074 

25 

PACS 100 

0.042 

0.018 

0.056 

25 

PACS 160 

0.028 

0.023 

0.066 

18 

SPIRE 250 

0.023 

0.015 

0.045 

15 

SPIRE 350 

0.017 

0.013 

0.036 

11 

SPIRE 500 

0.011 

0.007 

0.019 

6 



I. Bustamante et al.: Identification of new transitional disk candidates in Lupus with Herschel 


Fig. 2: SEDs of the 5 transitional disk candidates identified with Herschel. Blue dots represent the dereddened ancillary data from 
the literature from previous studies. The red diamonds are the clear detections with Herschel , and the red triangles the upper limits 
wherever no clear detection in the maps was found. The gray solid line is the median SED from the Lupus sample of the 29 non 
transitional disks in the region, and the gray shaded area the first and fourth quartiles. See Table[7]for more information. The dashed 
lines are the photospheric NextGen models from Allard et al. (2012). 
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I. Bustamante et al.: Identification of new transitional disk candidates in Lupus with Herschel 
Fig. 3: SEDs of the 29 YSOs with Herschel detections that do not fulfill our transitional disk selection criteria 
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Fig. 4: Fig. [3] continued. 
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I. Bustamante et al.: Identification of new transitional disk candidates in Lupus with Herschel 


Fig. 5: Images of each of the 34 sources with at least one point source detected by Herschel , plus object 2MASSJ1608537, added for 
completeness. The images correspond to a box of 120” x 120” in size. Coordinates are given in Table [5] The color scale is defined 
with the RMS of the pixel values (as background level) and this number plus three times the standard deviation (as maximum level). 
North is up and East is left. A circle indicates the expected position of the target. The standard deviation computed around the 
sources is shown in the bottom right side of the images. 
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Fig. 6: Fig. [5] continued. 
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I. Bustamante et al.: Identification of new transitional disk candidates in Lupus with Herschel 


Fig. 7: Fig. [5] continued. 
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I. Bustamante et al.: Identification of new transitional disk candidates in Lupus with Herschel 


Fig. 8: Fig. [5] continued. 






SSTc2dJ160901.4-392512 







~ 

■ .1 

1 . 

*■1 

100/im 

. i 

■ 

!■" 





* 

V ' J 


^ • 

-4 1 , 


V 

- 

V-T 

V 

( 


r 




V 

1 

* 

1 

■_ ; 

-■ 

y* p 

1 j 

u 

"" 

■■ - 

L- 

> 1 


Li 


rr: 0.038 


p rt. 

AN- 

\ --h-H — i 






500/ 

im 







r 





V 

y 






cr: 0. 

059] 




15 

























































































































































































































































































































































































































































































































































































I. Bustamante et al.: Identification of new transitional disk candidates in Lupus with Herschel 


Fig. 9: Fig. [5] continued. 
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